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Abstract 

The Bialas-Bzdak model of elastic proton-proton scattering assumes a 
purely imaginary forward scattering amplitude, which consequently van¬ 
ishes at the diffractive minima. We extended the model to arbitrarily large 
real parts in a way that constraints from unitarity are satisfied. The re¬ 
sulting model is able to describe elastic pp scattering not only at the lower 
ISR energies but also at yjs =7 TeV in a statistically acceptable manner, 
both in the diffractive cone and in the region of the first diffractive min¬ 
imum. The total cross-section as well as the differential cross-section of 
elastic proton-proton scattering is predicted for the future LHC energies 
of y/s =13, 14, 15 TeV and also to 28 TeV. A non-trivial, significantly non¬ 
exponential feature of the differential cross-section of elastic proton-proton 
scattering is analyzed and the excitation function of the non-exponential 
behavior is predicted. The excitation function of the shadow profiles is 
discussed and related to saturation at small impact parameters. 
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1 Introduction 


In a pair of recent papers the Bialas-Bzdak model (BB) [1] of small angle elastic 
proton-proton (jpp) scattering at high energies was studied at y/s = 7 TeV 
center-of-mass LHC energy [2,3]. 

In this work, we extend those investigations by improving on the original 
BB model by adding a real part to its forward scattering amplitude (FSA) in a 
unitary manner, and, furthermore, we present the extrapolation of the BB model 
to future LHC energies as well. Our method to include the energy evolution 
of the parameters is similar to the so-called “geometric scaling” discussed in 
Refs. [4] and [5], however, in our case only the constant and the linear terms 
are used as a function of In (s), while the quadratic terms used in Refs. [4] and 
[5] were not needed at our current level of precision. 

During 2014, the TOTEM experiment made public an important preliminary 
experimental observation at y/s = 8 TeV LHC energy: the pp elastic differential 
cross-section shows a deviation from the simplest, exponential behavior at low-f, 
where t is the squared four-momentum transfer of the pp scattering process [6]. 
This feature of the y/s = 8 TeV preliminary TOTEM dataset was related to 
t-channel unitarity of the FSA in Ref. [7], a concept that we also focus on, using 
and generalizing the quark-diquark model of Bialas and Bzdak to determine the 
shape of the FSA in elastic pp scattering. 

In its original form, the BB model [1] assumes that the real part of the 
FSA is negligible, correspondingly, the FSA vanishes at the diffractive minima. 
At the ISR energies of y/s =23.5—62.5 GeV, that were first analyzed in the 
inspiring paper of Bialas and Bzdak [1], this assumption is indeed reasonable, 
as it is confirmed in Ref. [2], At these ISR energies, only very few data points 
were available in the dip region around the first diffractive minimum of elastic pp 
scattering, which were then left out from the BB model fits of Ref. [2] to achieve 
a quality description of the remaining data points. However, in recent years, 
TOTEM data [9] explored the dip region at the LHC energy of y/s = 7 TeV in 
great details, at several different values of the squared four-momentum transfer 
t. Ref. [2] demonstrated, that the original BB model cannot describe this dip 
region, not without at least a small real part that has to be added to its FSA 
in a reasonable way. 

Subsequently, the BB model has been generalized in Ref. [3] by allowing for 
a perturbatively small real part of the FSA, which improved the agreement of 
the model with TOTEM data on elastic pp scattering at the LHC energy of 
y/s = 7 TeV. It was expected that the main reason for the appearance of this 
real part is that certain rare elastic scattering of the constituents of the protons 
may be non-collinear thus may lead to inelastic events even if the elementary 
interactions are elastic. The corresponding phenomenological generalization of 
the Bialas-Bzdak model [3] was based on the assumption that the real part of the 
FSA is small, and can be handled perturbatively. The resulting a-generalized 
Bialas-Bzdak (aBB) model was compared to ISR data in Ref. [3], and it was 
demonstrated that a small, of the order of 1 %o real part of the FSA indeed results 
in excellent fit qualities and a statistically acceptable description of the data in 
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the region of the diffractive minimum or dip. However, at the LHC energy 
of yfs = 7 TeV, although the real part of the fit becomes significantly larger 
than at ISR, the same aBB model does not result in a satisfactory, statistically 
acceptable fit quality, although the visual quality of the fitted curves improve 
significantly as compared to that of the original BB model [3]. 

These results indicate that at the LHC energies the real part of the FSA 
may reach significant values where unitarity constraints may already play an 
important role. The unitarity of the ^-matrix provides also the basis for the 
optical theorem, which in turn provides a method to determine the total cross- 
section from an extrapolation of the elastic scattering measurements to the t = 0 
point. In the aBB model [3], unitarity constraints were not explicitly considered: 
as the original BB model with zero real part obeyed unitarity, adding a small 
real part may result only in small (and expected) deviations from unitarity and 
the optical theorem. However, when the model was fitted to the 7 TeV TOTEM 
data in the dip region [3], the extrapolation to the point of t = 0 and the related 
value of the total cross-section underestimated the measured total cross-section 
by about 40%, suggesting, that perhaps the real part of the FSA may be large, 
and unitarity relations should be explicitly considered. 

These indications motivate the present manuscript, where the BB model 
is further generalized to arbitrarily large real parts of the FSA, derived from 
unitarity constraints. The resulting model is referred to as the real extended 
Bialas-Bzdak (ReBB) model. 

In Refs. [2] and [3] a terse overview is reported about the field of elastic 
scattering at high energies, that summarizes the developments up to 2013. In 
this introduction let us highlight only some more recent works, in order to put 
our results to a broader context of recent, independent investigations of elastic 
pp collisions at the LHC energies. 

The recent analysis of Ref. [16] applies a quark-diquark representation of 
the proton to describe TOTEM data [9] also using antiproton-proton data. The 
quark-diquark approach of Ref. [16] is a continuation of studies from the late 
sixties [14], which was first applied to describe pp elastic scattering in Ref. [15]. 
The early studies [14,15] provide the foundation of the quark-diquark represen¬ 
tation of the BB model as well. A recent improvement on this idea is Ref. [17], 
which introduces the idea of “Pomeron elasticity” to increase the real part of 
the FSA, which approach shows an interesting relationship with our present 
analysis. 

The structure of this manuscript is as follows: in Section 2, the general form 
of the FSA is re-derived for the case of a non-vanishing real part starting from 
^-matrix unitarity. Then this result is applied to the extension of the BB model 
to a non-vanishing and possibly large real part of the FSA. The mathematical 
relations between the resulting ReBB and the earlier aBB models are specified 
in subsection 2.2. 

In Section 3, the specified ReBB model is fitted to TOTEM data on elastic 
pp scattering at yfs = 7 TeV, both in the diffractive cone [10-12] and in the dip 
region [9]. 

Based on these fits and comparisons of the ReBB model to yfs = 7 TeV 
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data, the shadow profile function A(b) is evaluated in Section 4.1. This shadow 
profile function characterizes the probability of inelastic pp scattering at a given 
impact parameter b , and is compared to the shadow profile functions of elastic 
pp collisions at lower, ISR energies. Section 4.2 is devoted to study the structure 
of the differential cross-section da/dt at low-|t| values and also to compare it 
with a purely exponential behavior. 

In Section 5, the excitation function of the fit parameters is investigated and 
their evolution with yfs is obtained based on a geometrical picture. The model 
parameters are extrapolated to the expected future LHC energies of y/s =13, 14 
and 15 TeV, as well as for 28 TeV, that is not foreseen to be available at man¬ 
made accelerators in the near future, but may be relevant for the investigation of 
cosmic ray events. The excitation functions of the shadow profile functions A(b) 
are also discussed. Finally we summarize and conclude. 


2 The real extended Bialas-Bzdak model 

Although the original form of the Bialas-Bzdak model neglects the real part of 
the FSA in high energy elastic pp scattering, the model is based on Glauber 
scattering theory and obeys unitarity constraints. 

The phenomenological generalization of the Bialas-Bzdak model [3] is based 
on the assumption, that the real part of the FSA is small, and can be handled 
perturbatively. However, it turned out that the addition of a small real part 
does not lead to a statistically acceptable description of TOTEM data on elastic 
pp collisions at yfs = 7 TeV. In this manuscript, we consider the case, when the 
real part of the FSA is not perturbatively small. We restart from S-matrix 
unitarity, and consider how the BB model can be extended to significant, real 
values of the FSA while satisfying the constraints of unitarity. 

2.1 S-matrix unitarity in the context of elastic pp scatter¬ 
ing 

In this subsection some of the basic equations of quantum scattering theory 
are recapitulated. The scattering or S matrix describes how a physical system 
changes in a scattering process. The unitarity of the S matrix ensures that the 
sum of the probabilities of all possible outcomes of the scattering process is one. 
The unitarity of the scattering matrix S is expressed by the equation 

SS^ = I, (1) 

where I is the identity matrix. The decomposition S = I + iT , where T is the 
transition matrix, leads the unitarity relation Eq. (1) to 

T - T* = *TT t , (2) 

which can be rewritten in the impact parameter b representation as 

21m t e i(s,b) = |t eJ (s,6)| 2 + CT me z(s, b), (3) 
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where s is the squared total center-of-mass energy. 

The functions &i ne i{s,b) = d 2 ai ne i/d 2 b and \t e i(s, b)\ 2 = d 2 a e i/d 2 b are the 
inelastic and elastic scattering probabilities per unit area, respectively. The elas¬ 
tic amplitude t e i(s, b) is defined in the impact parameter space and corresponds 
to the fth partial wave amplitude Ti(s) through the relation £+1/2 by/ s/2, 
which is valid in the high energy limit, y/s —> oo. 

The unitarity relation (3) is a second order polynomial equation in terms of 
the (complex) elastic amplitude t e i{s, b). If one introduces the opacity or eikonal 
function [18-23] 


t e i(s,b) = i 


l _ e —f2(s,fe) 


di ne i can be expressed as 




Re f2(s,b) 


(4) 

(5) 


The formula for t e i is the so called eikonal form. From Eqs. (5) the real part of 
the opacity function f l(s, b) can be expressed as 

Re fi(s, b) = — i In [1 — d ine i (s, &)] . (6) 

In the original BB model it is assumed that the real part of t e i vanishes. In this 
case Eqs. (4) and (6) imply that 


t e i{s,b) = i 


l-V 7 ! 


&inel (- 5 , 6 ) 


(7) 


If the imaginary part Im Q is taken into account in Eq. (4) the result is 


tel(s,b)=i -ZinelM 


( 8 ) 


where the concrete parametrization of Imf2(s, b) is discussed later. 

To compare the model with data the amplitude Eq. (8) has to be transformed 
into momentum space 


T(s, A) 


+ OO +OO 

= I I e iA \i(s,b)d 2 b 

— OO —OO 


= 2wi / Jq (A • b ) 


l _ e -n(s,b) 


bdb, 


(9) 


where b = |6|, A = |A| is the transverse momentum and Jo is the zero order 
Bessel-function of the first kind. In the high energy limit, y/s —> 00 , A(f) ~ 
\f—t where t is the squared four-momentum transfer. Consequently the elastic 
differential cross-section can be evaluated as 

^ = T|T(,,A)|V (10) 
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According to the optical theorem the total cross-section is 

a t ot = 2 ImT(s, A)| t=0 , 

while the ratio of the real to the imaginary FSA is 

ReT(s, 0) 

P ~ ImT(s,0) ' 

2.2 The Bialas-Bzdak model with a unitarily extended 
amplitude 

The original BB model [1] describes the proton as a bound state of a quark 
and a diquark, where both constituents have to be understood as “dressed” 
objects that effectively include all possible virtual gluons and qq pairs. The 
quark and the diquark are characterized with their positions with respect to the 
proton’s center-of-mass using their transverse position vectors s q and Sd in the 
plane perpendicular to the proton’s incident momentum. Hence, the coordinate 
space H of the colliding protons is spanned by the vector h = (s q , s d , s q , s' d ) 
where the primed coordinates indicate the coordinates of the second proton. 

The inelastic pp scattering probability d’i ne i(s,b) in Eq. (7) is calculated as 
an average of “elementary” inelastic scattering probabilities a(h; b) over the 
coordinate space H [24] 


( 11 ) 

( 12 ) 


+oo +oo 


Zineiib) = (cr(h;b)^ ^ = J ... J d hp(h) ■ a{h;b), (13) 


— oo —oo 


where the weight function p{h) is a product of probability distributions 


P{h) = D(s q ,s d ) ■ D(s' q ,s' d ). (14) 

The D(s q , s d ) function is a two-dimensional Gaussian, which describes the center- 
of-mass distribution of the quark and diquark with respect to the center-of-mass 
of the proton 


D (Sg, s d ) 


1 + A 2 


o-{sl + < 


)/ R «“6 2 {s d + As 9 ), A 


m q 
m d ' 


(15) 


The parameter R qd , the standard deviation of the quark and diquark distance, 
is fitted to the data. Note that the two-dimensional Dirac S function preserves 
the proton’s center-of-mass and reduces the dimension of the integral in Eq. (13) 
from eight to four. 

Note that the original BB model is realized in two different ways: in one 
of the cases, the diquark structure is not resolved. This is referred to as the 
p = ( q , d) BB model. A more detailed variant is when the diquark is assumed to 
be a composition of two quarks, referred as the p = (g, (g, g)). Our earlier studies 
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using the aBB model indicated [3], that the p = ( q , d) case gives somewhat 
improved confidence levels as compared to the p = (g, (g, g)) case. So for the 
present manuscript we discuss results using the p = (g, d) scenario only, however, 
it is straightforward to extend the investigations to the p = (g, (g,g)) case. 
We have performed these calculations but we do not detail their results here 
given that they result in fits which are not acceptable at y/s = 7 TeV. To 
demonstrate that this p = (g, (g,g)) model does not work, we report only its 
rather unsatisfactory fit quality when the data analysis is discussed. 

It is assumed that the “elementary” inelastic scattering probability cr(h\ b) 
can be factorized in terms of binary collisions among the constituents with a 
Glauber expansion 


<r(h;b) = 1 - nil f 1 “ a o,b(b + s' a - s b ) , 

a b 


a, be {q,d}, 


(16) 


where the indices a and b can be either quark q or diquark d. 

The <j a b (s) functions describe the probability of binary inelastic collision 
between quarks and diquarks and are assumed to be Gaussian 

(Jab (s) = A ab e~ s2/s -\ S 2 ab = R 2 a + Rl, a,b G {q, d} , (17) 


where the R qi Rd and A a b parameters are fitted to the data. 

The inelastic cross-sections of quark, diquark scatterings can be calculated 
by integrating the probability distributions Eq. (17) as 


&ab, inel 


+oo +oo 

J j (Tab (s) cl 2 S = 7T AabSlb . 


—oo —oo 


(18) 


In order to reduce the number of free parameters, we assume that the diquarks 
are bounded very weakly, hence the ratios of the inelastic cross-sections a a b ,inel 
satisfy 

(Tqq, inel ■ (Tqd, inel ■ (Tdd : inel —1:2:4, (19) 

which means that in the BB model the diquark contains twice as many partons 
than the quark and also that these quarks and diquarks do not “shadow” each 
other during the scattering process. This assumption is not trivial. The p = 

(g, (g,g)) version of the BB model allows for different cr^mei : <r q d,inei ■ <Tdd, inel 
ratios. However, as it was mentioned before, the p = (g, (g, g)) is less favored 
by the data as compared to the p = (g, d) case presented below. 

Using the inelastic cross-sections Eq. (18) together with the assumption Eq. (19) 
the A q d and Add parameters can be expressed with A qq 

4 R 2 a 4 Rl , 

Aqd = A qq „ , Add = ' (^0) 

Kq + R d K d 

In this way only five parameters have to be fitted to the data R q d, R q , Rd, A, 
and Aqq. 
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The last step in the calculation is to perform the Gaussian integrals in the av¬ 
erage Eq. (13) to obtain a formula for di ne i(b). The Dirac 5 function in Eq. (15) 
expresses the protons’ diquark position vectors as a function of the quarks po¬ 
sition 

Sd — A Sq : Sd — A Sg . (21) 

After expanding the products in the Glauber expansion Eq. (16) the follow¬ 
ing sum of contributions is obtained 

(j(h, 6) —O qq T 2 • (Jqd A &dd i^&qq&qd A (Jqd A &qq&dd T 2 ‘&qd&dd') 

+ (, a qq IJ qd + 2 (J qq ,J qdC r dd A Vdd&qd) ~ ^qq^qd^dd > (22) 

where the arguments of the cr a b(s) functions are suppressed to abbreviate the 
notation. 

The average over H in Eq. (13) has to be calculated for each term in the 
above expansion Eq. (22). Take the last, most general, term and calculate the 
average; the remaining terms are simple consequences of it. The result is 

+oo +oo 

I = (-(Tqqa 2 q d (Jdd) H = / / dh P (K) ■ CT qq (T 2 qd CT d d) , (23) 


where the p{h ) weight function Eq. (14) is a product of the quark-diquark dis¬ 
tributions, given by Eq. (15). Substitute into this result Eq. (23) the definitions 
of the quark-diquark distributions Eq. (15) 



(24) 


where v = (1 + A 2 )/(2 • R 2 qd ) and the integral over the coordinate space H 
is explicitly written out; it is only four dimensional due to the two Dirac S 
functions in p(h). Using the definitions of the a a b (s) functions Eq. (17) and the 
expression A = A qq A q dAd q Add the integral Eq. (24) can be rewritten, to make 
all the Gaussian integrals explicit 


1 = - 


4v 2 A 


+oo +oo 


d Sqd~Sqe 


2 „/ „-2v( 


<+< 


^nib 


,-c k i(b-s k +s{) 2 


(25) 


—oo —oo 


where the abbreviations cu = Sk? refer to the coefficients in Eq. (17). Finally, 
the four Gaussian integrals have to be evaluated in our last expression Eq. (25), 
which leads to 


I 


L -b 2 S- 

-e B , 


B 


(26) 




where 


B — Cqd,dq (v T Cqq T ^ 6 dd ) T (1 Dqd,dq 



(27) 


and 



4i> + (1 + A)“ ( Cki + c mn ) , 

V ( Ckl + Cmn ) + (1 + A) CklC mn . 


(28) 


Each term in Eq. (13) can be obtained from the master formula Eq. (26), by 
setting one or more coefficients to zero, Cu = 0 and the corresponding amplitude 
to one, A k i = 1. 

Up to now, we have evaluated the real part of the opacity or eikonal function 
U(s, 6), that is determined by the inelastic scattering probability per unit area 
according to Eq. (6). Now we also have to specify the imaginary part of the 
complex opacity function, that determines the real part of the FSA. Here several 
model assumptions are possible, but from the analysis of the ISR data and the 
first studies of the 7 TeV TOTEM data at LHC we learned, that the real part of 
the FSA is perturbatively small at ISR energies, it becomes non-perturbative at 
LHC but the scattering is still dominated by the imaginary part of the scattering 
amplitude. 

We have studied several possible choices. One possibility is to introduce 
the imaginary part of the opacity function so that it is proportional to the 
probability of inelastic scatterings, which is known to be a decreasing function 
of the impact parameter b. A possible interpretation of this assumption may be 
that the inelastic collisions arising from non-collinear elastic collisions of quarks 
and diquarks follow the same spatial distributions as the inelastic collisions of 
the same constituents 


Im tt(s,b) = -a - dinei (s,b ), 


(29) 


where a is a real number, corresponding to the shape parameter of the differ¬ 
ential cross-section of elastic pp scattering. 1 

The above proportionality Eq. (29) between Imf2(s,&) and cq rae i(s,&) pro¬ 
vided the best fits from among the relations that we have tried but it is far from 
being a unique possibility for an ansatz. Note that the a = 0 case corresponds 
to the p = (q,d) version of the original BB model of Ref. [1], where the FSA 
has a vanishing real part. 

1 It is interesting to note a similarity of this concept with the so called “Pomeron elasticity”, 
which was introduced independently for quark-diquark models in Refs. [16] and [17]. Pomeron 
elasticity has a similar shape modifying role in the dip region: it was introduced to increase 
the real part of FSA by modifying the standard Pomeron trajectory in elastic quark-quark, 
quark-diquark and diquark-diquark scattering, that contributes significantly in particular close 
to the dip region of pp elastic scattering. 
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In the a < 0.1 perturbative limit the aBB model of Ref. [3] can also be 
obtained as follows. The aBB model can be defined by the relation 


Im 12(s, b) 


Ct * (Jinel b^) 
O'inel (-5, 6) 2 


(30) 


This definition is equivalent to the original form of the definition of the aBB 
model in Ref. [3] where <Ji ne i(s,b) in Eq. (7) was allowed to have a small 
imaginary part, in the form of cq ne j(s, b) —> (1 + ia)<ii ne i(s , b ). This condition 
was satisfied for fits at ISR energies, with a « 0.01, however, at LHC energies 
the aBB model did not describe the TOTEM data in a satisfactory manner [3]. 

We have also investigated the assumption that the real and the imaginary 
parts of the opacity function are proportional to one another 


Imfi(s, b) = — a ■ Ref2(s, b). (31) 

However, as the results using Eq. (31) were less favorable than the results 
obtained with Eq. (29), we do the data analysis part, described in the next 
section, using Eq. (29). 

We mention this possibility to highlight that here some phenomenological 
assumptions are necessary as the ReBB model does allow for a broad range of 
possibilities for the choice of the imaginary part of the opacity function. 

In this way, the ReBB model is fully defined, and at a given colliding energy 
only six parameters determine the differential (10) and total cross-sections (11) 
and also the p parameter, defined with Eq. (12). The parameters that have to 
be fitted to the data include the three scale parameters, R q , Rd , R q d, that fix 
the geometry of the pp collisions, as well as the three additional parameters a, 
A and A qq . Two of the latter three can be fixed: A = 0.5 if the diquark is very 
weakly bound, so that its mass is twice as large as that of the valence quark. In 
practice we also fix A qq = 1, assuming that head-on qq collisions are inelastic 
with a probability of 1, according to Eq. (17). 

Thus in the present data analysis two out of six parameters are fixed and only 
four parameters are fitted to the data at each i/s: the three scale parameters 
R qi Rd and R q d , as well as the shape parameter a. As we shall see, the shape 
parameter a will play a key role when describing the shape of the dip of the 
differential cross-sections of elastic pp scatterings at LHC energies. 


3 Fit method and results 

The pp elastic differential cross-section data measured by the TOTEM exper¬ 
iment at 7 TeV is a compilation of two subsequent measurements [9-11]. The 
squared four-momentum transfer value t sep = —0.375 GeV 2 separates the two 
data sets. 2 Note, that the two datasets were taken with two very different 
settings of the machine optics of the LHC accelerator. 

2 The squared four-momentum transfer value t 3ep separates the bin centers at the common 
boundary, the two bins actually overlap [9—11]. 
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The ReBB model, defined with Eqs. ( 8 ) and (29), was fitted to the data 
at ISR energies and at LHC energy of ^fs = 7 TeV. The relation between the 
imaginary part of S2(s, b) and a is defined with Eq. (29). In agreement with our 
previous investigations the A qq = 1 and A = ; parameters can be kept constant, 
which reduces the number of free parameters to four R q d, R q , Rd and a. The 
case when the parameters A qq and A are also free is summarized at the end of 
this section. 

In the course of the minimization of the ReBB model we take into account 
the uncertainty of the overall scale factor of the measured data. The fitted y 2 
function to the data points is 



(32) 


where N is the number of fitted data points, dcq/dt is the ith measured dif¬ 
ferential cross-section data point and da t h,i/dt is the corresponding theoretical 
value at the ith data point calculated from the ReBB model. The value cq is the 
mere statistical uncertainty of the ith data point. We set ai um i to a conservative 
value of 5 % at each ISR energy[27]. In case of the TOTEM data sets it is set 
to 4 %, which is determined by the uncertainty of the CMS luminosity[10,11]. 
Parameter 7 is the additional normalization parameter to be minimized and 
70 = 1 stands for the value of luminosity scaled to unity, while ai um i is the 
relative luminosity uncertainty reported by the measurement. 

The sum in Eq. (32) runs over the fitted data points, and the additional 
term takes into account the contribution of the luminosity uncertainty [25]. 

First we attempted to fit the ReBB model in the 0 < |t| < 2.5 GeV 2 range, 
fitting simultaneously both the low-|f| TOTEM measurement of Ref. [11] and 
the one containing the dip region[10], see Fig. 1. Note, that in this particular 
fit two normalization parameter were used: 71 below and 72 above |t aep |, since 
the mentioned two TOTEM measurements are independent. This fit provides 
X 2 /NDF = 289.04/158 = 1.83 and CL = 9 x lO " 9 < 0.1%, which is statis¬ 
tically not an acceptably good fit quality, although, as indicated in Fig. 1, the 
fit looks reasonably good by eye. So, unfortunately, we could not get a unified 
and statistically acceptable description of the differential cross-section of elastic 
pp scattering in the whole measured t interval in the framework of the unitarily 
extended Bialas-Bzdak model. 

It is important to note, that we determined the fit quality using the statis¬ 
tical and the luminosity uncertainty only, where the latter is a t-independent 
systematic uncertainty. According to the original TOTEM publications [10,11], 
the systematic uncertainties of the two TOTEM data sets are very different due 
to different data taking conditions, especially due to the different LHC machine 
optics. The f-dependent part of the systematic errors allow the data points, as a 
function of |f|, to be slightly moved in a correlated and t dependent way, which 
could, in principle, improve our fit quality. 

However, this part of the systematics is rather difficult to handle correctly in 
the present analysis. So we decided to analyze the two TOTEM data sets sepa- 
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rately. If the separated fits to y/s = 7 TeV elastic differential cross-section dcr/dt 
data are evaluated, below and above the separation \t sep \, quality results can 
be obtained, which are shown in Figs. 2 and 3, and are summarized in Table 1. 
As detailed below, this strategy leads to reasonable fit qualities (CL = 1.8 %, 
statistically acceptable fit in the cone region and CL = 0.04 %, statistically 
marginal fit in the dip region), with a remarkable stability of fit parameters. 3 


y/s [GeV] 

23.5 

30.7 

52.8 

62.5 

7000 

1*1 [GeV 2 ] 

(0,2.5) 

(0, |tsep|) 

(|t„ P |,2.5) 

X^/NDF 

124.7/101 

95.6/46 

96.1/47 

76.2/46 

109.9/81 

120.4/73 

CL [%] 

5.5 

2 x 10~ 3 

3 x 10~ a 

0.3 

1.8 

4 x 10 -2 

R q [ fm\ 

0.27T0.01 

0.28T0.01 

0.28T0.01 

0.28±0.01 

0.45±0.01 

0.43T0.01 

Rd [ fm\ 

0.72T0.01 

0.74T0.01 

0.74T0.01 

0.75±0.01 

0.94T0.01 

0.91T0.01 

R q d [. fm ] 

0.30±0.01 

0.29T0.01 

0.31T0.01 

0.32±0.01 

0.32±0.05 

0.37T0.02 

a 

0.03±0.01 

0.02T0.01 

0.04±0.01 

0.04±0.01 

0.11±0.04 

0.12T0.01 

7 

1.01±0.05 

0.98T0.05 

0.90±0.06 

0.97T0.05 

1.00±0.04 

1.00 (fixed) 


Table 1: The values of the fitted ReBB model parameters from ISR to LHC 
energies. At the 7 TeV LHC energy, the pp elastic dcr/dt data measured by the 
TOTEM experiment is a composition of two subsequent measurements, which 
are separated at t sep = —0.375 GeV 2 . The errors and the values are rounded 
up to two valuable decimal digits. 

It is quite remarkable, that although the low-|f| fit ends at t sep = —0.375 
and the dip position is at significantly larger values of |t|, still the fit when 
extrapolated to the dip and the Orear region after the dip reproduces the data 
very well [13]. 

The calculated total cross-section of the low-|£| fit is crtot — 99.3 ± 3.8 mb, 
where the uncertainty is the propagated uncertainty of the fit parameters, in¬ 
cluding that of the shape parameter that can be relatively badly determined 
in the cone region. This result nevertheless agrees very well with the value 
Gtot = 98.0 ± 2.5 mb measured by the TOTEM experiment at y/s = 7 TeV [12] 
in a luminosity independent way, and the calculated value for the parameter 
p = 0.09 ± 0.03 is also reasonable, as within its errors it is consistent with the 
measured value of p = 0.145±0.091 as reported by the TOTEM experiment [12]. 
Note, that the uncertainty of the value of our p parameter is the propagated un¬ 
certainty of the ReBB model fit parameters as it is given in Fig. 2, consequently, 
it contains the effects of propagated statistical and luminosity uncertainties only. 

On the other hand, if we look at the fit to the dip region, with 7 = 1 
fixed, we also see a remarkable stability of the shape of the differential cross- 
section at low values of |f| that still yield reasonable values for the total cross- 
section (fj tot = 91.9 ±2.6 mb) and similarly reasonable values for the parameter 

3 The aBB version of the ReBB model [3] has already been fitted to ISR data but only 
in the restricted 0.36 < |t| < 2.5 GeV 2 range in order to be consistent with the /.-range of 
the available TOTEM data set of that time [10], while in this work, the low |t| data are also 
included at each energies. The more limited fit range explains the seemingly better y 2 /A T DF 
and CL values reported in our earlier publication Ref. [3] using the crBB model. 
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p = 0.10±0.01. The stability and consistency of the model description is visible 
also in Figs. 2 and 3. 

For the sake of completeness, we present also one of our fits at the ISR energy 
of y/s = 23.5 GeV, as indicated in Fig. 4 [9,11,26,27]. The parameters of the 
best fits and the parameters’ errors at each analyzed ISR energy are summarized 
in Table 1. 
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p+p -> p+p, diquark as a single entity at \fs=7000.0 GeV 



Figure 1: The fit of the ReBB model in the whole 0 < |t| < 2.5 GeV 2 range 
at y/s = 7 TeV. The fit uses the statistical errors of the data points and the 
luminosity error of the systematic uncertainty according to Eq. (32). Although 
the fit quality is not satisfactory, CL -C 0.1%, the fit looks good by eye. The 
fitted parameters are shown in the left bottom corner, parameters without errors 
were fixed during the MINUIT optimization. The total cross-section atot and the 
parameter p are derived quantities according to Eqs. (11) and (12), respectively. 
Parameter values are rounded up to two decimal digits. The uncertainty of the 
MINUIT error matrix after the fit is 2.8 %. 
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p+p -» p+p, diquark as a single entity at \fs=7000.0 GeV 



Figure 2: The fit of the ReBB model in the 0 < \t\ < \t sep \ range at y/s = 7 TeV. 
The fit quality is satisfactory, CL > 0.1%. The fit uses the statistical errors of 
the data points and the luminosity error of the systematic uncertainty according 
to Eq. (32). The fitted curve is shown with solid line, its extrapolation above 
\t se p | is indicated with a dashed line. The extrapolated curve remains close to 
the data points, following the measured differential cross-sections well even far 
away from the fitted region. The fitted parameters are shown in the left bottom 
corner, rounded to two valuable decimal digits. The parameters without errors 
were fixed in the minimization. The total cross-section atot and the parameter 
p are derived quantities according to Eqs. (11) and (12), respectively. The 
uncertainty of the MINUIT error matrix after the fit is 1.5 %. 
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p+p -» p+p, diquark as a single entity at \fs=7000.0 GeV 



Figure 3: The same as Fig. 2, but the fit is performed in the \t sep \ < |t| < 
2.5 GeV 2 range. The fit uses the statistical errors of the data points and we 
present the results for the 7 = 1 fixed case. The fitted curve is shown with solid 
line, its extrapolation is indicated with a dashed line. Parameter values are 
rounded up to two valuable decimal digits. Note that when the curve is extrap¬ 
olated to the low-|t| region, the extrapolated curve again follows the measured 
differential cross-section remarkably well even far away the fit region: the ReBB 
model fit is remarkably stable over the whole |t|-range. 


16 













p+p -> p+p, diquark as a single entity at \fs=23.5 GeV 



Figure 4: The fit of the ReBB model at t/s = 23.5 GeV in the 0 < |f| < 2.5 GeV 2 
squared four-momentum transfer |f| range. The fit uses the statistical errors of 
the data points and the luminosity error of the systematic uncertainty according 
to Eq. (32). Parameter values are rounded up to two valuable decimal digits. 

If the two parameters A qq and A are released and included into the set of 
fitted parameters of the ReBB model, the fit quality improves at each analyzed 
energy from the point of view of mathematical statistics. In case of y/s = 
30.7 GeV the improvement is quite significant as the confidence level of the 
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fit reaches CL = 8% instead of 2 x 10 -3 %. However, the two new fitted 
parameters introduce more correlations, which lead to large fit uncertainties, 
and the parameters A and A qq within errors remain in the range of their values 
that were fixed for the fits reported in Table 1, however, due to the correlations 
between the fit parameters when A qq and A are released, the fit parameters 
fluctuate more when evaluated as a function of y/s, consequently their trend 
is more difficult to determine. Therefore, in order to determine the excitation 
function of the model parameters, we utilized the results of the ReBB fit results 
as listed in Table 1, where these two parameters A = 1/2 and A qq = 1 are fixed. 

Also note that if ImS2(s,6) is defined to be proportional to Ref2(s, 6), ac¬ 
cording to Eq. (31), the MINUIT fit result of y 2 /NDF = 504.9/159 = 3.2 is 
obtained at y/s = 7 TeV in the 0 < |t| < 2.5 GeV 2 range, which is disfavored as 
compared to fits with Eq. (29), see also Fig. 1. 

In our introduction we shortly mentioned the p = ( q , (q,q)) version of the 
RcBB model, when the diquark is assumed to be a composition of two quarks [3]. 
At y/s = 7 TeV in the 0 < |f| < 2.5 GeV 2 range this scenario provides a fit 
result with y 2 /NDF = 15509/159 ~ 97.5, which means that the p = (q, ( q,q )) 
ReBB version can be rejected. The failure of this version is basically due the 
wrong shape of the differential cross-section: the second diffractive minimum 
appears too close to the first one. 

4 Discussion 

4.1 Shadow profile functions and saturation 

The fits, from which the model parameters were determined, also permit us to 
evaluate the shadow profile function 

A(s, b) = 1 — |exp [— fl(s, 6)]| 2 . (33) 

The obtained curves to A(b) are shown in Fig. 5. The shadow profile functions 
at ISR energies exhibit a Gaussian like shape, which smoothly change with the 
center-of-mass energy yfs. In this case, the A(s, b = 0) < 1 value indicates that 
the protons are not completely “black” at the ISR energies, even at their centre 
they do not scatter with the maximum possible probability per unit area. At 
the LHC energy of yfs = 7 TeV something new appears: the innermost part of 
the distribution shows a saturation, which means that around 6 = 0 the shadow 
profile function becomes almost flat and stays close to A(b) ss 1. Consequently, 
the shape of the shadow profile function A(b) becomes non-Gaussian and some¬ 
what “distorted” with respect to the shapes found at ISR. At the same time, 
the width of the edge of the shadow profile function A(b), which can be visu¬ 
alized as the proton’s “skin-width”, remains approximately independent of the 
center-of-mass energy yfs. 
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Figure 5: The shadow profile functions A(b) indicate a saturation effect at LHC, 
while at ISR energies a Gaussian shape can be observed. Note that the 7 TeV 
(black dashed) curve is based on the statistically acceptable fit result in the 
0 < |t| < 0.38 GeV 2 range. The distributions’ edge shows approximately the 
same width at each energy, corresponding to a constant “skin-width” of the 
proton. 


4.2 Non-exponential behavior of da e i/dt 

To compare the obtained low-|t| ReBB fit of Fig. 2 with a purely exponential 
distribution the following exponential parametrization is used 

dcr e i _ dcr e i 
dt dt 

where dcx e i/dt| t=0 = 506.4 mb/GeV 2 and the slope parameter of B = 19.89 GeV -2 
is applied, according to the TOTEM paper of Ref. [11]. 

The result, shown in Fig. 6, indicates a clear non-exponential behavior of 
the elastic differential cross-section in the 0.0 < |t| < 0.2 GeV 2 range at yfs =7 
TeV. 

A change of the slope parameter B(t , s) = (d/dt) ln(d<r/dt) around —t ~ 0.10 
GeV 2 was reported already in the year 1972 in elastic pp collisions at the ISR 
energy range of 21.3 < y/s < 54 GeV and the 0.02 GeV 2 < — t < 0.40 GeV 2 
squared four-momentum transfer range. A very similar structure, a deviation 
from an exponential behavior was also reported as early as in 1984 in the analysis 
of proton-antiproton elastic scattering at 546 GeV by Glauber and Velasco [28]. 
The first preliminary results on such a non-exponential behavior at the CERN 
LHC energy of 8 TeV were made public recently in pp elastic scattering by the 
TOTEM experiment at various conference presentations during 2014 [6] and 
were already interpreted in the theoretical work of Ref. [7], as the consequence 
of two-pion exchange and f-channel unitarity, while Ref. [8], related this t de¬ 
pendence of the slope parameter B to pion loop contribution to the Pomeron 
trajectory in two-channel eikonal models. 
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Figure 6: The ReBB model fit result, shown in Fig. 2, with respect to the 
exponential fit of Eq. (34). In the plot only the 0.0 < |t| < 0.2 GeV 2 range 
is shown, but the ReBB model is fitted to 7 TeV TOTEM data in the low-t 
interval of (0,t sep ). The ReBB fit result indicates a significant deviation from 
the simple exponential at low-|t| values. 


5 Extrapolation to future LHC energies and be¬ 
yond 

The ReBB model can be extrapolated to energies which have not been measured 
yet at LHC. The fit results of Table 1 and the parametrization 

P(s)=po+Pi-la.(s/so) (35) 

is applied for each parameter P £ {R q ,Rd,R q d,&}, where so = 1 GeV 2 . The 
parametrization Eq. (35) implies that the four free parameters of the original 
ReBB model that were fitted at each colliding energy independently, corre¬ 
sponding to altogether 20 free fit parameters at the 5 energies analyzed in this 
manuscript are now replaced with eight parameters pi that prescribe their en¬ 
ergy dependence. These fits to the energy dependence of the ReBB model 
parameters are shown in Fig. 7 and the fit parameters are collected in Table 2. 

The logarithmic dependence of the geometric parameters on the center-of- 
mass energy y/s in the parametrization Eq. (35) is motivated by the so-called 
“geometric picture “ based on a series of studies [29-34], In case of the a pa¬ 
rameter, which is not a geometrical property of the proton, the logarithmic y/s 
dependence is an additional assumption. As indicated by Fig. 7 and Table 2, 
such an energy dependence of the a shape parameter is consistent with the 
currently available data. 

Table 2 shows that the rate of increase with y/s, parameter p \, is an order of 
magnitude larger for R q and Rd than for R q d. The saturation effect, described 
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Figure 7: The excitation function of the parameters of the ReBB model, col¬ 
lected in Table 1, is determined from fits with Eq. (35) to each of the parameters 
R qi Rd, R q d and a. The plots about the resulting fits are collected here, while 
the parameters of the excitation functions are collected in Tabic 2. The statis¬ 
tically acceptable quality of these fits allow the ReBB model to be extrapolated 
to center-of-mass energies which have not been measured yet at LHC. 


in Section 4.1, is consistent with this observation as the increasing components 
of the proton, the quark and the diquark, are confined into a volume which is 
increasing more slowly. 


Parameter 

R q [ fm} 

Rd [fm] 

Rqd [fm] 

a 

X 2 /NDF 

6.2/3 

2A/3 

7.5/3 

1.2/3 

CL [%] 

10.2 

49.4 

5.8 

75.3 

Po 

0.15 ±0.01 

0.59 ±0.01 

0.30 ±0.01 

-0.04 ±0.01 

Pi 

0.017 ±0.001 

0.019 ±0.001 

0.002 ±0.001 

0.009 ±0.001 


Table 2: The fits of the ReBB model to data at each ISR energy and at the 7 TeV 
LHC energy, summarized in Table 1, allow for an extrapolation of the model 
parameters as a function of the center-of-mass energy yfs. The parametrization 
Eq. (35) is applied to each parameter of the ReBB model and the fits are shown 
in Fig. 7. Numerical values are rounded up to two valuable decimal digits. The 
fit quality information is provided in the first and second row of the table. Note 
that the fit quality is acceptable for each parameter as CL>0.1 % for each fit. 
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Using the extrapolation formula 
Eq. (35) and the value of the pa¬ 
rameters from Table 2 it is straight¬ 
forward to calculate the values of 
the parameters at 8 TeV, where the 
TOTEM measurement of the total 
cross-section is published [35], and 
at expected future LHC energies of 
y/s =13, 14, 15 TeV and also at 28 
TeV, which is beyond the LHC capa¬ 
bilities. Using the extrapolated val¬ 
ues of the parameters we plot our 
predicted pp elastic differential cross- 
section curves at each mentioned en¬ 
ergy in Fig. 8. The shadow profile 
functions A(b) can be also extrapo¬ 
lated, see Fig. 9. The shadow pro¬ 
file functions even allow us to visual¬ 
ize the increasing effective interaction 
radius of the proton in the impact pa¬ 
rameter space in Fig. 10. 

It is also important to see how 
the most important features change 
with center-of-mass energy y/s: the 
extrapolated values of the total cross- 
section at.ot, the position of the first 
diffractive minimum \tdi P \ and the pa¬ 
rameter p is given in Table 3. 

Our calculated value at yfs = 8 TeV is atot = 99.6 mb, which is consistent 
with the total cross-section cr tot = 101.7 ± 2.9 mb at yfs = 8 TeV, measured 
with a luminosity-independent method by the TOTEM experiment [35]. 

According to Table 3, the predicted value of \tdi P \ and atot moves more than 
10% when yfs increases from 8 TeV to 28 TeV, while the value of C = \tdi P \- a to t 
changes only about 2 %, which is an approximately constant value, within the 
errors of the extrapolation. 



-t [GeV 2 ] 

Figure 8: The pp elastic differential 
cross-section is extrapolated to 8 TeV 
as well as to future LHC energies and 
beyond. 
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Figure 9: The shadow profile function at the extrapolated energies yfs. The 
results show the increase of the proton interaction radius with increasing yfs 
energies. Also note that the “edge” of the distributions remains of approxi¬ 
mately constant width and shape. 


[TeV] 

a to t [mb] 

tdzp [GeV^] 

P 

\tdi P \ ■ a tot [mb GeV" 1 ] 

8 

99.6 

0.494 

0.103 

49.20 

13 

106.4 

0.465 

0.108 

49.48 

14 

107.5 

0.461 

0.108 

49.56 

15 

108.5 

0.457 

0.109 

49.58 

28 

117.7 

0.426 

0.114 

50.14 


Table 3: The extrapolated values of the total cross-section a tot at future LHC 
energies and beyond. The position of the first diffractive minimum \tdi P \, the 
parameter p and the \tdi P \ ■ <Jtot value is also provided at each energy. Note that 
the predicted value of \tdi P \ and atot moves more than 10% when yfs increases 
from 8 TeV to 28 TeV, while the value of |idj P | • a to t changes only about 2 %. 


A similar, and exact, scaling can be derived for the case of photon scattering 
on a black disk, where the elastic differential cross-section is [36] 


deblock 

dt 


7 tR 4 


Ji(q-R) 

q ■ R 


(36) 


where R is the radius of the black disk and t = —q 2 . The total cross-section is 
given by 


^tot,black — 27 lR . (37) 

In this simple theoretical model the position of the first diffractive minimum, 
following from Eq. (36), and the total cross-section Eq. (37) satisfies 

Cblack — \tdi P ,black \ ' &tot 7 black — 27rJi,l(^m) ~ 35.9 mb GeV , (38) 
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where j' 1,1 is the first root of the first order Bessel-function of the first kind J\ ( x ). 

The scaling behavior, indicated by the stability of the value C, is observed, 
but it is significantly different from the black disk model, described by Eq. (38), 
as the corresponding value Cuack is significantly different 


Cblack 7 ^ C . 


(39) 


In this sense the value of C indicates a more complex scattering phenomena, 
than the scattering of a photon on a black disc, however, the constancy of the 
product suggests the validity of an asymptotic geometric picture, in agreement 
with the recent observations in Refs. [3,37-39]. 


b x [fm] 



Vs~=23 GeV 
b x [fm] 



by 


[fm] 


y[7 =14 TeV 


b x [fm] 



[fm] 


-3-2-10 12 3 


V7 =7 TeV 
b x [fm] 



b y [fm] 


-3-2-10 12 3 
V7 =28 TeV 


Figure 10: Visualization of the shadow profile functions A(b) in the transverse 
plane of the impact parameter vector (b x ,b y ). The figures show the increase 
of the proton effective interaction radius in the impact parameter space with 
increasing center-of-mass energy yfs. It can be also observed that the black 
innermost core of the distributions is increasing, while the thickness of the pro¬ 
ton’s “skin”, the gray transition part of the distributions, remains approximately 
independent of the center-of-mass energy yfs. 
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6 Summary and conclusions 

The real part of the forward scattering amplitude (FSA) is derived from unitarity 
constraints in the Bialas-Bzdak model leading to the so-called ReBB model. The 
added real part of the FSA significantly improves the model ability to describe 
the data at the first diffractive minimum. In total the ReBB model describes 
both the ISR and LHC data in the 0 < |t| < 2.5 GeV 2 squared four-momentum 
transfer range in a statistically acceptable manner; in the latter case the fit 
range has to be divided to two parts, according to the compilation of the two 
independent TOTEM measurements. The results are collected in Table 1. 

The fit results also permit us to evaluate the shadow profile functions A(b), 
see Fig. 5. The plots indicate a Gaussian shape at ISR energies, while at LHC 
a saturation effect can be observed: the innermost part of the shadow profile 
function A(b) around b = 0 is almost flat and close to A(b) ps 1. The elastic 
differential cross-section can be compared to a purely exponential distribution 
and the comparison shows a significant deviation from pure exponential in the 
0.0 < |t| < 0.2 GeV 2 range. 

The fit results allow the determination of the excitation functions of the 
ReBB model at future LHC energies and beyond, with parameters collected in 
Table 2 and predicted differential cross-section curves shown in Fig. 8. The 
shadow profile functions can be also extrapolated, see Fig. 9, which predicts 
that the saturated part of the proton is expected to increase with increasing 
center-of-mass energy y/s. The edge of the distribution, the “skin-width” of the 
proton, expected to remain approximately constant. It is worth to mention that 
the extrapolated version of the ReBB model utilizes of only eight parameters, 
the pi parameters of Table 2, and in this sense a “minimal “ set of parameters 
is applied. 
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